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Abstract 

This paper gives an introduction to and overview of the functional quantum programming language QML. 
The syntax of this language is defined and explained, along with a new QML definition of the quantum 
teleport algorithm. The categorical operational semantics of QML is also briefly introduced, in the form 
of annotated quantum circuits. This definition leads to a denotational semantics, given in terms of su- 
peroperators. Finally, an implementation in Haskell of the semantics for QML is presented as a compiler. 
The compiler takes QML programs as input, which are parsed into a Haskell datatype. The output from 
the compiler is either a quantum circuit (operational), an isometry (pure denotational) or a superoperator 
(impure denotational). Orthogonality judgements and problems with coproducts in QML are also discussed. 

Keywords: QML, language, functional, teleport, denotational semantics, operational semantics, Haskell. 



1 Introduction and motivation 

Language development for quantum computation is a rapidly developing research 
area [4] , motivated by the application of established formal reasoning and verifica- 
tion techniques within a quantum framework, understanding the behaviour of quan- 
tum computation, aiding the development of new algorithms and gaining a deeper 
understanding of how they work. This paper discusses the syntax and features of, 
and gives a compiler for, a language allowing both classical and quantum control: 
QML [1,5]. The syntax and semantics for QML is a complete redevelopment of that 
presented previously [1], as the language has been changed to remove a problematic 
interpretation of coproducts (section 2.2); the interpretation of orthogonality has 
also been updated (section 2.1). In addition, in this work the operational semantics 
is made concrete by a compiler for QML programs implemented in Haskell. 

QML has a syntax similar to other functional languages [4], and is based on 
strict linear logic: i.e. linear logic with contraction, but without implicit weakening 
(in contrast, Selinger's QPL uses affine linear logic [8], which allows weakening but 
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not contraction). QML also integrates reversible and irreversible computations in 
a single language, where weakenings (which can give rise to the collapse of super- 
positions and entanglement) must be explicit. Contractions are allowed and are 
modelled as a form of sharing, analogous to the behaviour of classical functional 
languages. Differences between QML and other languages include the use of a quan- 
tum control operation, provided by the quantum-if construct if°. To use the if° 
operation, the branches of the computation must be orthogonal (distinguishable), 
the proof of which is supplied automatically by the type checker at compilation, 
and hence the programmer need not supply the condition, nor need it appear in the 
syntax of terms. Classical control is provided by a second classical-if construct if, 
which measures the term being branched over. 

QML has both an operational and denotational semantics [5], supporting formal 
reasoning principles with an algebra for equational reasoning and a normal form [2] . 
The operational semantics of QML is presented using a categorical formalisation of 
the quantum circuit model, which is realised by a compiler that translates QML 
programs into 'typed' quantum circuits. The denotational semantics is also imple- 
mented, which translates QML terms, via the operational semantics category FQC, 
into either isometries Q° or superoperators Q. 

Recent developments in QML include a fully operational compiler, with a type 
inference algorithm for QML terms, and the automatic derivation and extension of 
the orthogonality judgements and circuits required for quantum control, which have 
all been implemented. The implementation of the orthogonality judgements is such 
that it can be easily expanded as new rules for proving the orthogonality of terms are 
added to QML. This short paper outlines the new syntax of QML and its semantics. 
A new and faithful interpretation of the quantum teleportation algorithm in QML 
is also presented as an example of using QML. Details for using the compiler, and 
the output it generates, are then provided. Finally, the future development of the 
language, the semantics and compiler, and the uses of the compiler are described. 

2 The syntax of QML 

This section introduces the syntax for QML (see also ref. [5]). The symbols a, r 
are used to vary over QML types, given by a = Q\ \ Q2 \ cr ® r, where the type 
constructor is the tensor product <g) corresponding to a product type and Q2 is a 
qubit type, x, y, z are used to vary over names. Typing contexts (T, A) are given 
by r = • j r, x : a where • is the empty context. Contexts correspond to functions 
from a finite set of variables to types. 

Constants k,leC are also used to define the syntax of expressions. Function 
variables are used to refer to previously defined QML programs. The terms of QML 
consist of those of a first-order functional language, extended with quantum data, 
a quantum control structure, and a measurement operator. The vector notation y 
is used for sequence variables to be measured (weakened). The grammar of QML 
terms is defined thus: 

(Variables) x,y,... £ Vars 

(Prob amplitudes) k, l, ... 6 C 

(Patterns) p,q ::= x \(x,y) 
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( Terms) t, u ::= x \ x v | () | (t, u) | let p = t in u 

if t then u else u' | if° i then u else m' 
qfalse y | qtrue y | k x t \ t + m 
Quantum data is modelled using the constructs k x t and i + u. The term k x t, 
associates the probability amplitude n with the term t. The term t + u describes a 
quantum superposition of t and u. Quantum superpositions are first class values, 
and can be used in conditional expressions to provide quantum control. For example: 
if ° (qtrue + qfalse) then t else u evaluates both t and u and combines the results 
in a quantum superposition. Note that the term Ao x t + Ai x u, where t, u are not 
qubits, is syntactic-sugar for if° (Ao x qtrue + Ai x qfalse) then t else u. The 
type-checker and orthogonality judgements ensure that this is a valid operation, by 
providing a proof that t and u are orthogonal (distinguishable in some way), that 
their types match, and that they are strict terms (they produce no garbage). 

In a quantum control operation, the two branches must be orthogonal, other- 
wise the type system would accept terms that implicitly perform measurements. 
Without this restriction "valid" programs could be written in QML that are not 
physically realisable by a quantum computer. Orthogonality judgements are in- 
ferred automatically by static analysis of terms (see section 2.1). 

As an example of superposition formation, the term (-i=) x qfalse + (-^) x qtrue 
is an equal superposition of qfalse and qtrue. Normalisation factors that are equal 
may be omitted. 

Finally, a QML program is a sequence of function definitions, where a function 
definition is given by / V = t : r. A Haskell-style syntax is used to present program 
examples. For example, the QML function below (left) is equivalent to the following 
Haskell- like code (right): 

/ (xi : ai,x 2 :a 2 ,-.,x n :a n ) f : o\ -o a 2 . . -o a n — o t 
= t : t f xi x 2 . . x n = t 

2. 1 QML orthogonality judgements 

QML has a basic type system that tracks the use of variables, preventing them from 
being weakened inappropriately. However, the type system still accepts terms which 
implicitly perform measurements. As a consequence QML would accept programs 
which are not realisable as quantum computations. 

Consider the expression if° x then qtrue else qtrue. This expression returns 
qtrue without using any information about x. In order to maintain the invariant 
that all measurements are explicit, the type system should reject the above expres- 
sion. More precisely, the expression if° x then t else u should only be accepted 
if t -L u. This notion intuitively ensures that the conditional operator does not 
implicitly discard any information about x during the evaluation. The branches of 
a superposition should also be orthogonal for similar reasons. 

Mathematically, two terms, t, u, are orthogonal if their inner-product is equal to 
zero, (t\u) = 0. If this is the case then the judgement t _L u is true, but if the inner- 
product yields any other value then t is not orthogonal to u. In the presentation of 
an equational theory for QML [2] the orthogonality judgements are replaced by an 
inner-product judgement on terms, to much the same effect. However, the inner- 
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product approach is more informative and flexible, and gives a method of reasoning 
about orthogonality, hence in future this method may be adopted for all terms. 
The following rules give the current QML orthogonality judgements: 



qtrue _L qfalse qfalse X qtrue 

t _L u t _L u 

pair _L paii^ 



(t, v) _L («, w) (v, t) _L (w, u) 

t _L u Aq«o = — A^ki 
Ao x i + Ai x bIko x i + Ri x ii 
t _L u Aq^o = —XiK>i 



sup 



if° (Aq X qtrue + Ai X qfalse) then t else u _L if° (kq x qtrue + Ki X qfalse) then t else u 

t _L u t _L u t ± u t ± u' 



supif° 



t _L if° c then u else u' if° c then u else «' _L t 



if? 



The first two axioms state that the basic states of qtrue and qfalse are orthogonal. 
The third and fourth rules state that pairs of terms can be orthogonal, provided that 
one component of a pair is orthogonal to the other pair's corresponding component. 
The two sup rules state when superpositions of terms can be orthogonal; the second 
is a restatement of the first, translating superpositions using if°. The final two rules 
state that if° statements can be orthogonal if all the component terms are. 

The use of if ° in QML programs is valid only if the two branches are orthogonal. 
Hence, for the Hadamard operation (section 4), it is required that — qtrue+qfalse _L 
qtrue + qfalse, with the appropriate renormalisation. In this case the _L sup rule 
verifies orthogonality. 

The rules for orthogonality given so far are incomplete, and may be extended. 
Orthogonality judgements must also be interpreted by the operational semantics, 
discussed in ref [5]. 



2.2 Removing coproducts from QML 

In a previous version of QML [1], here referred to as QML®, the language included 
the notion of a tensorial coproduct, denoted by ©. This coproduct has now been 
removed. The types of QML® were generated by Qi, a<S)T, and <r©r. Qubits were 
not primitive, but defined as Q2 = Qi ffi Si- The coproduct allowed any finite type 
to be directly represented in this way; not just limited to 0,2- The introduction 
rules used for © were the usual coproduct rules of a left and a right injection: 

T\- a s:cr r\- a t:r 

hintroi Mntro r 

r h a inl s : a ® r r h a inr t : cr r 

The coproduct type was interpreted ascrffir = Q2©l cr U T li where \a U r| could 
store a value of either \a\ or |r|, by padding the smaller type. Using the coproduct 
and injection rules, qfalse and qtrue were defined in QML® as fnl() : Q2 and 
qfalse = inr() : Q2, respectively, omitting the weakening property of QML®. 

Instead of if and if° rules, QML® implemented two ©-elimination rules: case, 
providing classical-control (a generalisation of if), and a quantum-control operation 



4 



Grattage 



case (generalising if ). The quantum (non-measuring) ©-elimination rule is similar 
to the standard coproduct elimination rule, and is given as: 

r \- a c : o-0r 
A, x : a h° t : p 



r®Ah° case c of {inl x => t | inr y => u} : p 

The non-strict (measuring) case removes the orthogonality requirement, and 
does not require sub-terms to be strict. The if a rules (a £ {o, — }; if a = o then 
the rules are strict, i.e. measurement-free) would then be derived as: 

if a b then t else u = case a b of {inl _ =^ t | inr _ =4> u} 
The branches of a case operation can be of different sizes, and this was dealt 
with in the semantics of QML® by padding the type of the smaller branch. The 
padding of one type in this way can lead to the garbage becoming entangled with the 
useful output in some way. This happens, for example, if branching over Q\ (g) Q2. 
The garbage created by padding may indirectly measure the qubit which is being 
branched over. Consequently, this approach is not compositional, and has therefore 
been rejected. 

This version of QML resolves this issue by removing the coproduct. Qubits 
are now primitive, as are if and if°. Additionally, the strict if° does not allow 
any garbage to be produced. Coproducts may be reintroduced, possibly limited to 
classical types. 



3 An operational semantics for QML 



The new operational semantics for QML is briefly discussed here, which is an up- 
dated version of that presented in [1]. This semantics is defined by giving a trans- 
lation from QML terms into morphisms in the category FQC. FQC morphisms 
consist of a reversible quantum circuit (f), the input context T, the output type a, 
and the size of the auxiliary heap h and any garbage g. Any heap qubits are initially 
set to {false) in the computational basis, and garbage qubits can be removed by 
the partial trace operation at the end of the computation. A full development of 
FQC is given in references [5,6]. 

As an example, the classical if construct is defined by the following typing rule 
and operational semantic: 



r h c : Q 2 
Aht.u : a 



■ if 



c e fqc r q 2 

t e FQC A a 
u S FQC A a 



A h if c then t else u 



where 



iF 0p c t u e FQC (r (g) A) a 
lF 0p c t u = Oc + h c + h t \ u ,g c 

in the operational semantics is the following circuit: 



<t>t\4>u 



PC 



with (pc as a "context-splitting" operation that copies any variables used by both 
subterms. 4>t\<Pu denotes a conditional circuit which performs cj>t if the result of the 
conditional circuit c is true (|1)) and (j) u if it is false (|0)). 
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In this way, the semantics of QML is defined recursively over the syntax of terms, 
and results in a valid FQC morphism for a valid QML term. Given a QML term, 
the QML compiler follows the operational semantics and outputs an FQC morphism 
represented as a typed circuit, which can then either be displayed, exported for use 
with other programs, or used to directly implement the program on a quantum- 
circuit based quantum computer. The FQC morphism can be further evaluated 
via the denotational semantics, producing a unitary matrix Q- (for strict-QML 
terms without heap), an isometry Q° (for strict-QML terms), or a superoperator Q 
(for QML terms that produce garbage). In addition, orthogonality judgements and 
circuits for quantum control and superpositions of terms are automatically inferred 
by the compiler at compile-time, so there can be no orthogonality errors during run- 
time. The relationships between QML without measurement (QML°), full QML 
(QML), typed quantum circuits (FQC), and isometries (Q°) and superoperators 
(Q), is shown in the following diagram: 




Q° r a> : *- Q r a 

where [-]op denotes the operational semantics, and [-]d denotes the additional 
translation from quantum circuits into the denotational semantics. The full seman- 
tics can be found in reference [5]. 

4 Example: Teleportation in QML 

The quantum teleportation describes how to transport a quantum state using a small 

amount of classical communication. A QML interpretation of the teleportation 

circuit with deferred measurement has previously been presented, along with a full 

description of the algorithm [5]. However, this circuit relies on the existence of a 

quantum channel. In order to demonstrate and explain the syntax of QML and the 

compiler, a new, faithful, implementation of the quantum teleport algorithm, which 

explicitly makes use of measurement, is presented. 2 This algorithm is similar to 

that developed by Selinger and Valiron [9], and also in reference [3], which includes 

a relevant discussion of teleportation, both with and without deferred measurement. 

These examples show the elegance of allowing quantum control via the quantum if° 

construct {CNot) and term superpositions (Epr), in writing functions. For simple 

one qubit functions (Had), the branches of the quantum control are the columns 

from the unitary matrices that describe the operation. A simple measurement 

operation, using classical control, is also defined in the following example which 

implements the teleportation algorithm (Tele): 
Had, Qnot, Meas e Q2 — o Q2 

Had b = if ° b then qfalse + — qtrue — The Hadamard operation 
else qfalse + qtrue 

Qnot b = if b then qfalse else qtrue — The Not operation (Pauli X) 

Meas b = if b then qtrue else qfalse — A measurement operator, using if 

CNot £ Q2 — o Q2 —o Q2 ® Q2 —A quantum cnot operation 

2 The code is included with the compiler as teleport . qml, see section 5. 
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CNot s t = if° s then (qtrue, Qnot t) 
else (qfalse, t) 

Epr&Q2®Q2 - The EPR pair, |00> + |11> 

Epr = (qtrue. qtrue) + (qfalse, qfalse) 

Bmeas S Q2 —o Q2 — o Q2 C*D Q2 — The Bell-measurement operation 

Bmeas x y = let [x 1 ', j/') = CNot x y 

in (Meas (Had x'), Meas y') 
U S Q2 — o ® Q2 — o QS — The unitary correction operations 

U q xy = let (x, 2/) = sj/ in if x then (if y then Un q else Uio q) 

else (if y then f/ni </ else g) 

C/oi,C/io,C/n &Q2 ^Q2 
Uqi x = if ° x then qfalse else qtrue 
U\o x = if ° a; then — qtrue else qfalse 
Un x = if° x then — qfalse else qtrue 

— The quantum teleportation algorithm 
Tele e Q2 Q2 

Tele q = let (a, b) = Epr — a is given to Alice, b to Bob 

/ = Bmeas q a — Result of Alice's Bell-measurement is classical data 

in U b f — Bob applies U to his qubit, using classical data / 



5 The QML compiler 

This section briefly describes the design and operation of the Haskell QML compiler, 
which can be found on the project website. 3 The objective is to take a file containing 
a QML program, consisting of QML function definitions, and output an annotated 
(typed) quantum circuit which realises the QML program as an FQC object. This 
is achieved by implementing the operational semantics in Haskell. Additionally, 
the circuit produced by the compiler can be further processed to produce either a 
unitary matrix (Q~), isometry (Q°), or superoperator (Q), as appropriate, giving 
an implementation of the denotational semantics of QML, as factored through the 
operational semantics. 

The compiler has a modular design, giving a clear logical structure. For example, 
the QOrth module contains all the code for generating the orthogonality judgements 
and circuits, while QCirc contains the definition of the circuit datatypes and asso- 
ciated functions. The compiler exploits advanced Haskell features, such as monads 
and pattern matching, and making use of the ideas put forward by Vizzotto et 
al [10]. The operational semantics is realised in the QComp (QML Compilation) 
module. 

To compile QML functions into the operational semantics (FQC morphisms, 
represented as typed circuits), the "run typed circuit" command is used: 
runTC "#f ilename" "#f unction". For example, to evaluate the typed circuit of 
the function Epr (from section 4) the command is runTC "teleport . qml" "Epr". 
This outputs the following typed circuit (as a Haskell datatype): 

^-{H\— f — 2 2 
l \X}- 2 2 

where there is no input context, the heap is two qubits (marked h), and a pair 



3 Instructions for downloading the QML compiler can be found on the QML compiler website 
http://fop.cs.nott.ac.uk/qml/compiler. The Haskell compiler GHC is also required. To use the com- 
piler, the QML system must be loaded into GHC via the command "ghci qml". This will initialise the 
system with all the modules required to interpret and evaluate QML programs. 
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of qubits are the output. H and X are the Hadamard and Pauli-X operations. 
As heap is initialised to qfalse, this circuit describes the function 1 00) — > 1 00) + 
1 11), producing an EPR pair. The compilation of a QML term into a circuit is 
efficient; each QML term is recursively translated directly into an FQC morphism 
(an annotated circuit), as described in section 3. No quantum computation is 
simulated, so this does not effect the efficiency of the compiler - it is a simple 
recursive translation into circuits. The output is further optimised after compilation, 
by collapsing identities and permutations and other simple circuit manipulations 
which do not effect the action of the circuit (see QCirc). 

There are three main options for further evaluation of a QML term: 

• The QML term could be evaluated to a unitary matrix (Q~), interpreting only the 
reversible circuit <j) from the FQC morphism. As this option classically computes 
the full mathematical interpretation of the program it is inherently inefficient. 
The output is actually a triple (h, g, cj>) £ (Int, 7ni,Unit), where h, g gives the size 
of any heap required or garbage produced. The command for producing a unitary 
matrix is runM; 

• The function could be evaluated to an isometry (Q°), which initialises any re- 
quired heap, and is the full description for terms that produce no garbage. This 
option is no less efficient than the runM option, and is the preferred evaluation 
option. The output is actually a pair (g, cp) € (Int, Q°), where g gives the size of 
any garbage (if the QML function is impure) . The command for evaluating to an 
isometry is runl; 

• Thirdly, the QML term can be interpreted as a superoperator (Q), which ini- 
tialises heap and traces out any garbage, using the command runS. This option 
is substantially less efficient than the previous two options, as the state space is 
doubled and the partial trace is a computationally expensive operation. 

Together, the options runl and runS give an interpretation of the denotational 
semantics of QML factored through the category FQC, as shown in the diagram 
in section 3. A direct implementation of the denotational semantics, without using 
the operational semantics, is an extension currently being developed. Please refer 
to the project website for full details. 

6 Conclusions and further work 

This paper introduces the language QML and presents its semantics with a compiler 
and example programs. The semantics and compiler give a realisable interpretation 
of QML programs as quantum circuits in a formal, categorical, setting. The seman- 
tics can be extended in many ways, such as expanding the current orthogonality 
judgements, or by the addition of non-linear, classical, data. The algebra for the 
pure fragment of QML [2] is currently being extended to include measurement, fol- 
lowing work on van Tonder's quantum lambda calculus [3]. It would be instructive 
to implement this algebra, especially the normal form, as part of the QML compiler. 
Future possibilities for the development of the language also including developing a 
notion of higher-order functions for QML, and adding iteration to the language. 
The development of QML and the compiler is an ongoing project which has 
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already reached a functional state. As the language and semantics evolve, exten- 
sions and new features can be incorporated into the compiler; which also provides 
a useful testbed for the development of new language features and capabilities. For 
example, an extension of the orthogonality circuits given in [5] was developed using 
the compiler in this way. The compiler also facilitates the testing and development 
of new QML algorithms, such as the described teleportation algorithm. It has also 
been useful in allowing others to experiment with quantum programming and get 
immediate feedback on the behaviour of their functions, in a style that is famil- 
iar to computer scientists, logicians, and physicists with functional programming 
experience. 

Further extensions to the compiler include adding the ability to export typed 
circuits as images, or in notation compatible with tools such as MatLab and Math- 
ematica. Possible relationships with the measurement calculus, the Haskell QIO 
monad [7], and other formalisms are being studied, and may provide new insights. 
This will lead to new features being developed, such as basis independence, and 
further useful abstractions. 
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